Scattering theory of nonlinear thermoelectric transport 
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We investigate nonlinear transport properties of quantum conductors in response to both electrical and thermal 
driving forces. Within scattering approach, we determine the nonequilibrium screening potential of a generic 
mesoscopic system and find that its response is dictated by particle and entropic injectivities which describe 
the charge and entropy transfer during transport. We illustrate our model analyzing the voltage and thermal 
rectification of a resonant tunneling barrier. Importantly, we discuss interaction induced contributions to the 
thermopower in the presence of large temperature differences. 
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may dramatically enhance S in strongly correlated systems as 
in magnetically diluted metallic hosts 



Introduction. Recent advances in nanoscale thermoelec- 
tric materials suggest novel functionalities and highly im- 
proved performances [1]. A key ingredient of thermoelectric 
devices is the Seebeck effect, which depends on the simulta- 
neous existence of thermal and electric driving forces. As a 
result, energy conversion from waste heat is possible under 
the conditions of zero net current. The Seebeck coefficient 
S measures the amount of thermovoltage generated across a 
conducting sample when a thermal gradient is externally ap- 
plied. Interestingly, the thermoelectric figure-of-merit is pro- 
portional to S 2 . Therefore, it is highly desirable to put for- 
ward new routes to increase S. Electron-electron interactions 

H and artificial Kondo 

impurities 

On the other hand, large temperature drops give rise, quite 
generally, to thermal rectification effects [4]. The possibility 
to apply sharp thermal gradients seems to be more feasible in 
nanostructured materias, as recently demonstrated in superlat- 
tices with periods spanning a few nanometers |5t]. Strikingly 
enough, a self-consistent theory of nonlinear thermoelectric 
transport valid for quantum conductors is still lacking. This is 
the gap we want to fill in this work. 

Linear thermoelectric effects within the scattering approach 
were discussed in Ref. 6. At the same time, pioneering ex- 
periments analyzed the main properties of the thermopower 
at linear response in quantum point contacts [7] and quantum 
dots Q8|] . Subsequent advances have unveiled fluctuating ther- 
mopower in chaotic dots [9], large S in Andreev interferome- 
ters II 1 Ofl and thermoelectric anisotropies in multiterminal bal- 
listic microjunctions [11]. The Seebeck coefficient can also 
help determine the conduction character of a molecular junc- 
tion [12]. Only recently has been possible a clear observation 
of thermal rectification effects in mesoscopic systems 11311 . 
Thus, it is natural to ask how phase-coherent current and ther- 
mopower are affected in the nonlinear regime of transport. 

In the isothermal case, all terminals are held at the same 
background temperature T. Refs. 03 and [l5| then provide a 
convenient theoretical framework to include nonequilibrium 
effects beyond linear response. The theory is based on an 
expansion around the equilibrium point but, importantly, the 
nonlinear transport coefficients are complicated functions of 



the screening response of the conductor out of equilibrium. 
This purely interaction driven response is described in terms 
of characteristic potentials which measure how the internal 
potential counterbalance the ensuing charge pile-up due to 
a voltage shift. Hence, the characteristic potentials depend 
on the particle injectivity of those carriers originated in the 
shifted terminal. The role of these particle injectivities is 
crucial because they determine departures from the Onsager- 
Casimir symmetry relations 116, 17] ubiquitously found in 
nonlinear transport experiments I 1 8h22I1 . Here, we show that 
when the system is perturbed with a temperature shift its re- 
sponse is dictated by entropic injectivities which quantify the 
entropy transported in the charge imbalance process. Below, 
we discuss the role of both particle and entropic injectivities in 
two conceptually simple but generic problems — the formation 
of rectifying terms in thermally driven electric currents and 
the differential Seebeck coefficient beyond linear response. 

Theoretical model. We consider a mesoscopic conduc- 
tor coupled to multiple terminals a,/3... characterized with 
bias voltages eV a = \i a — Ep is the electrochemical 
potential and Ep the Fermi energy) and temperature shifts 
Oa = T a — T (T a is the reservoir temperature). The elec- 
tronic transport is completely determined by the scattering 
matrix s a p = s a p[E,eU{r)\ which, in general, is a func- 
tion of the carrier energy E and the potential landscape in- 
side the conductor U(r) ffl [l5tl . In turn, U(f) is a func- 
tion of position r and the set of voltage and temperature 
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s a pS a p], the electrical cur- 



shifts. Defining A a p 

rent is expressed as I a = ^ J dEA a p(E)fp{E) where 
fp(E) = 1/(1 + exp [{E - Ep - Vp)/k B Tp\) is the Fermi 
distribution function in reservoir /?. In the weakly nonlinear 
regime of transport, the dominant terms appear up to second 
order in an expansion of the electric current in powers of the 
driving fields V a and 6 a : 
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The electrical and thermoelectric linear conductances are 

H G a p = -(2e 2 /h) J dEA aP d E f * (2e 2 /h)A aP {E F ) 
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and L aP = -(2e/hT)JdE(E - E F )A aP d E f =s 
(eir k%T/3h) dEA a p\E=E F > respectively, where the approx- 
imate expressions correspond to a Sommerfeld expansion to 
leading order in fc^T '/Ep. Here, / is the Fermi distribution 
function when all V a and 6 a are set to zero. We empha- 
size that the linear conductances are evaluated at equilibrium 
and, as a consequence, G a p and L a p are independent of the 
screening potential U. The situation is completely different 
for the nonlinear coefficients. We find, 
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where E a p = [(E — E F )/T]d E A a p. Notably, the nonlin- 
ear responses depend on how the scattering matrix changes, 
through the potential U, in response to a shift in voltage or 
temperature. Since we are concerned with small changes 
away from equilibrium, an expansion of U up to first order 
suffices: 



(3) 



where u a = (dU/dV a )eq and z(dU/d9 a ) eq are characteris- 
tic potentials that describe the internal change of the system 
to a shift of voltage and temperature, respectively, applied to 
terminal a. In the sequel, we derive the self-consistent proce- 
dure to determine the electrostatic potential in the presence of 
electrical and thermal forces. 

The net charge response of the system away from its equi- 
librium state can be decomposed into two terms, namely, the 
bare charge injected from lead a and the screening charge that 
builds up in the conductor due to interaction with the injected 
charges: q = (/bare + <7scr- The contribution to qbare due to a 
voltage imbalance in lead a is given by the particle injectiv- 
ity i/P(E). This is a partial density of states associated with 
scattering states that describe those carriers originated from 
lead a [14]. In addition, a shift of temperature in lead a also 
induces a change in q\, are . In contrast to the voltage case, how- 
ever, where every carrier with an energy E contributes pos- 
itively to gbare> in the thermally bias case the contribution of 
a temperature shift in lead a gives rise to a heat addition or 
removal depending on whether the carrier energy E is larger 
or smaller than Ep 112311 . This crucial fact must be reflected in 
the entropic injectivity denoted by v^: 



*f3a 



dE 



(4) 



E 



ds 



- E F | 
T Sf3a ~dE 



0a 



(5) 



To be concise, we have assumed that the potential is homo- 
geneous (i.e., position-independent) within the sample (the 
extension to inhomogeneous fields is straightforward 01511 ) 
and that the WKB approximation applies in order to make 
the replacement 6/SU — > —ed/dE. We note that the fac- 
tor (E — E F )/T represents the entropy transfer associated to 
adding a single carrier 112411 . Then the accumulation or de- 
pletion bare charge imbalance due to voltage or to tempera- 
ture shifts becomes qtare = eE Q K e ^ + D e a Qa) where 
DP = - J dEuP(E)d E f, and^ - - f dEis e a (E)d E f. 
Next, we obtain the screening charge from the response of 
the internal potential, AU = U — C/ eq , to changes in the leads' 
chemical potential and temperature. Within the random phase 
approximation, one has q sa — e 2 HAU. II is the Lindhard 
function which in the static case (frequency-dependent effects 
are not considered here) and in the long wavelength limit reads 
II = - J2 a D a = -D at T = [D = D(E F ) is the sample 
density of states] [25]. These approximations are excellent for 
our purpose since (i) if T ^ one can simply replace the pre- 
vious expression with II = J dED(E)d F f and (ii) the long 
wavelength limit amounts to carrier energies well below the 
tunnel barrier heights that couple the conductor to the exter- 
nal reservoirs. But this is precisely the range of validity of the 
WKB approximation used to express D p and D e in terms of 
energy derivatives only. 

Our set of equations is closed when we relate the out-of- 
equilibrium net charge with AU employing the Poisson equa- 
tion, V 2 AU = -Airq. We use Eq. (f3]l and the fact that V a and 
8 a shifts are independent. We then identify a pair of separated 
equations: 



-V 2 U a 

-V 2 z„ 



4ne 2 Uu a = A^e 2 D p a 
- 4ire 2 ILz a = A-KeD e a . 



(6) 
(7) 



These equations become nonlocal in the case of inhomoge- 
neous fields. 

The voltage and temperature derivatives, dg A a p and 
dv A a p, can be determined once the characteristic potentials 
are known since dg A a p = z 1 8A a p/5U — > —ez^dpAap 



and d V -,A a f3 = u 7 SA al3 /SU 
becomes 



y d F A a p. Thus, Eq. (ffl 
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d E f, 
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This is our central result. Importantly, Eq. ([8]) is not only of 
formal interest but offers clearly practical advantages. 

Quantum dot. As an illustrative application of the formal- 
ism exposed above, we now investigate the nonlinear thermo- 
electric transport properties of a quantum dot when Coulomb 
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FIG. 1. (Color online) Self-consistent screening potential U for a 
quantum dot system with Vi = E F + V/2, V2 = E F - V/2, T\ = 
T+6/2, and T 2 = T-0/2 (see inset). Parameters: T = 0.05, Ti = 
T2 =0.1 (we take ks = 1). (a) U at 6 = and various dot level 
positions: Ed — — 1 (black solid line), Ed = (green dashed line), 
and Ed — — 1 (blue dot line). The red dot-dashed line corresponds 
to the leading order approximation U = uV = r/V/2 independently 
of the dot level position, (b) U at V = and E d = ±1, 0. Colored 
dotted lines correspond to U = z9 with the characteristic potential 
calculated from Eq. d 14b . 



interactions are treated within a mean-field approximation. 
Preliminary observations suggest interesting nonlinear ther- 
moelectric effects in quantum dots [26]. We consider a single 
level with energy Ed coupled to two reservoirs (1 and 2) via 
tunnel barriers (see inset of Fig.[TJ. Thus, the level adquires a 
broadening given by T = Ti + The corresponding Breit- 
Wigner line shape depends, quite generally, on the internal po- 
tential U, which is to be self-consistently calculated through 
the Poisson equation. The dot charge is then 



qd -*J dE (E-E d -eUy+Ti 



(9) 



We expand Eq. (O to leading order in V a , a and U, and 
subtract the equilibrium dot charge q e d [given by Eq. (O with 
./': h /I: 

5q d = e 2 D p V 1 + e 2 D P V 2 + eD\Q x + eD e 2 9 2 - e 2 DU , (10) 

with Sq d = qd — q d denoting the charge excess due to voltage 
and temperature shifts and 
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T (E- E d ) 2 + r 2 
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d E f- (12) 



Clearly, Eqs. ( fTTT i and ( fl2b agree with the integrated parti- 
cle and entropic injectivities of Eqs. © and © if the Breit- 
Wigner representation of the scattering matrix is used. 

In a discrete form, the Poisson equation is written in terms 
of a geometrical capacitance C which connects electrically 
the dot to an external gate terminal. Accordingly, the charge 



excess of the dot obeys Sqg = C(U — V g ) where V g is the gate 
potential. Using Eq. ( fTOb we find the dot internal potential 



U = 



e 2 D p 2 V 2 



eDf 0i + eD e 2 6 2 + CV g 



C 



(13) 



from which the characteristic potentials follow, 
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1(2) 



C + e 2 D 



(14) 

Rectification effects. We consider the charge neutral limit 
(C = 0) since it applies to the experimentally relevant case 
of strong interactions. Moreover, if the dot is symmetrically 
biased (Vi = E F + V/2, V 2 = E F - V/2, 7\ = T + 6/2, and 
T 2 = T - 0/2) then u = dU/dV = r//2 and z = dU/dO = 
(Df - Dl)/[2e{D p + D p )] to leading order in V and with 
r/ = (Ti — T 2 ) /T the tunneling asymmetry lfl6ll . In Fig.[TJwe 
show the exact dot potential obtained from a numerical, self- 
consistent calculation of Eq. (O compared to its approximate 
value [Eq. (U3H, We distinguish between the isothermal case 
[0 = 0, Fig. [T](a)] and the electrically unbiased case [V = 0, 
Fig. [T](b)]. In tne former, the self-consistent potential is plot- 
ted for three values of the dot level Ed — ±1, 0. The curves 
for the fully self-consistent potential U agree with approxima- 
tion U = uV at low V only, as expected. At higher voltages 
in the strongly nonlinear regime and for Ed = ±1, higher or- 
der terms (V 2 or higher) become relevant and U departs from 
its linear dependence. We recall, however, that linear response 
transport coefficients depend on U cq only and they are insen- 
sitive to the variation of U with V. Only the nonlinear current 
allows us to explore this regime. Interestingly, at resonance 
(Ed = 0) the contributions to U from even powers in V are 
absent. In the electrically unbiased case [see Fig. [TJb)], we 
present U in response to a thermal shift for Ed = ±1,0. Par- 
ticularly interesting is the particle-hole symmetry case E d = 
for which the internal potential vanishes to all powers. We 
also compare the full calculation with the leading-order ap- 
proximation U = z6. Notice that contrary to the electrically 
biased case z depends on the dot level position. The agree- 
ment is quite reasonable at low temperature shifts, although at 
larger higher order terms become comparable. 

The evolution of the current for an electrically and ther- 
mally driven quantum dot is shown in Fig. |2|a) and Fig. |2jb), 
respectively, for fixed Ed = 1. For = the current first 
follows Ohm's law at low V and then, at higher voltages, ac- 
quires a quadratic dependence in V; the latter produces rectifi- 
cation effects. The I—V curves can be approximated up to V 2 
with I = G\\V + GmV 2 + 0(V 3 ) where the leading-order 
nonlinearity in the Sommerfeld approximation reads 



Cm = -rd E An\E=E F (l - 2ui) . 
h 



(15) 



and depends on the internal potential response. The I—V 
curves in Fig. |2ja) correspond to three values of the tunneling 
asymmetry (black solid line for rj = 0, red dashed line for 
rj = 0.5, and blue dashed-dotted line for 77 = 0.95) and show 
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FIG. 2. (Color online) Electrical current for a quantum dot system 
with Vi = E F + V/2, V 2 = E F - V/2, 6»i = T -6/2, and 
02 — T + 6/2 for three different values of the tunneling asymmetry 
rj. (a) I-V characteristics for 9 — along with the leading-order 
nonlinearity / ~ GuV + GmV 2 . The latter correspond to the 
colored dotted lines, (b) 1-9 characteristics for V = 0. Colored 
dotted lines correspond to / ~ L\\9 + Lm# 2 . Upper inset: Lm 
versus T as a function of 77 for Ed = 1. Lower inset: Lm for 
E d = ±1,0 and 77 = 0. 

good agreement with the second order expansion (colored dot- 
ted lines) except for very high voltages where deviations due 
to higher order terms (V 3 and so on) occur. In Fig. |2jb) we 
show / driven by a temperature shift for V = 0. We compare 
the full 1-9 characteristics (for different 77 values) with the 
second-order expansion in 9, I = L\\8 + L\\\8 2 + 0(9 3 ), 
where the thermal rectification term is 

eTT 2 k 2 T 

Lm = — — {dEAn\E=E F — 2ziTd E An\E=E F ) , 

(16) 

to leading order in the Sommerfeld expansion. First, the cur- 
rent grows linearly with 9 and then higher orders in 9 become 
relevant as 9 increases above a threshold where Lm is large 
enough. We plot Lm in the upper inset of Fig. |2|b) and 
find a nonmonotonous behavior with the background temper- 
ature T. For completeness, we also show in the lower inset of 
Fig- Etb) the dependence of the Lm for various level posi- 
tions of a symmetrically coupled dot (77 = 0). Interestingly, 
in the particle-hole symmetry point both L\\ and Lm vanish 
identically whereas for E4 = ±1, Lm presents an opposite 
behavior as a function of T. 

Thermopower. The thermopower S yields the voltage 
generated across the sample in response to an applied ther- 
mal bias at vanishing current condition. In the linear trans- 
port regime and for a two-terminal conductor, the Seebeck 
coefficient is So = V/0\i=o = —L\\/G\\. This ex- 
pression is correct in the limit 9 — > 0. At low temper- 
atures, it can be approximated to the Mott formula So ~ 
(■7T 2 k 2 3 T/6)dE In A\\\e=e f oc T whereas for high T we find 
So ^ (-k B /e)(E d - E F )/T oc T~ l . In Fig. Ha) we nu- 
merically calculate So for an electrically biased quantum dot 
(Vi = — V% = —V/2) when only one reservoir is heated 
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FIG. 3. (Color online) (a) Linear-response thermopower So for a 
symmetrically voltage biased dot and one heated contact (9i — 9 and 
#2 = 0) at Ed = 1. Low and high temperature limits of So are ex- 
plicitly shown, (b) Thermopower S beyond linear response for three 
different background temperature values. We show with colored dot- 
ted lines the leading-order expansion S — So + Si9 calculated from 
the sensitivity given by Eq. dl7t . 



(9i = 9, and 62 = 0). So has the correct L-dependence 
at low and high temperature. More interesing are the 9- 
corrections to S when 9 is not small. Then, we can expand 
S = S + Si9 + 0(9 2 ) where the Si is the thermopower 
sensitivity which measures the deviations of S from a con- 
stant value. Importantly, a measurement of the differential 
thermopower dS/ d9 gives precisely Si to leading order in 9. 
Specializing Eq. (fT} to the two-terminal case and setting 1 = 
we find 

Si = —^-[GuiL^+LuiG^ + GnLuiM^-Mui)], 

(17) 

valid when a single lead is heated. Inserting Eq. (ffj) in Eq. ( TFTb 
we compare the sensitivity with an exact calculation of S for a 
quantum dot as above. We observe in Fig. [3a) that excellent 
agreement is found for low 8 and that departures depend on 
the particular value of T. 

Conclusions. We have presented a general nonlinear scat- 
tering theory for mesoscopic conductors that are driven by 
electrical and thermal gradients. In the weakly nonlinear 
regime, screening effects arise in response to charge pile-up 
due to voltage or temperature differences. Importantly, the 
transmission probability becomes a function of the thermal 
gradient. We have found that the screening response can 
be described in terms of particle and entropic injectivities. 
We have illustrated our theory with an application to a two- 
terminal quantum dot setup, evaluating the current-voltage 
and current-temperature characteristics. Importantly, we have 
discussed thermopower sensitivity in the nonlinear regime of 
transport. Our result is relevant in view of recent advances in 
thermoelectrics at the nanoscale. 
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